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Abstract 

The Maslov correction to the wave function is the jump of — 7r/2 in the phase when the 
system passes through a caustic point. This phenomenon is related to the second variation 
and to the geometry of paths, as conveniently explained in Feynman's path integral frame- 
work. The results can be extended to any system using the semiclassical approximation. 
The 1-dimensional harmonic oscillator is used to illustrate the different derivations reviewed 
here. 

arXiv : quant -ph/0702236 

1 Introduction 

A fascinating fact, first observed at the end of the 19th century [1 , is that, after passing through 
a focal point, the phase of light jumps by —71/2. If a light beam is split into two parts one of 
them passes through a focal point while the other does not, then, when the two partial waves 
are recombined, a destructive interference is observed. 

This curious phenomenon has later been extended to the quantum mechanics of massive 
particles, where it is referred to as the "Maslov correction" [11 |3l II]- A simple illustration is 
provided by the harmonic oscillator [5] . Feynman's path integral framework [6j is ideally suited 
to understand how this comes about [7j. 

The aim of this Review is to derive the Maslov correction and study some if its aspects in 
the path integral framework. 

2 The Feynman propagator of the oscillator 

Let us first consider, for simplicity, a problem in 1 space dimension. Let xi and X2 two points 
and let T > a time interval]^ Let us indeed consider the set, denoted by V of all curves 7(t) 
between xi and 2:2, i. e. such that 7(0) = xi and 7(T) = X2- The classical action is a real 
valued function defined on P, 




(2.1) 
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^It is tacitly assumed that the system is conservative. Were this not the case, one should consider two instants 
ii and t2 > ti instead of T — t2 — ti. 



where L(x, x) is the Lagrange function of the system. Then the Principle of Least Action 
(Hamilton's Principle) tells us that the actual motion which starts in xi and arrives to X2 in T 
is the 7(i) that makes S stationary. 

Apart of exceptional cases, (see below) the two given points xi to X2 are joined, in time T, 
by a unique classical motion curve 7(t). But in Quantum Mechanics the situation is different. 
Intuitively, a quantum particle moves not along a unique path, but along all paths which join 
these points in the given time ! Feynman, in his Thesis, argued indeed that to any path 7 in 7-" 
is associated a complex number, namely 



exp 



(2.2) 



where S denotes the classical action (2.1) calculated along 7. Next, if the "amplitude" (whose 
square \ip\'^ is the probability) of finding our particle in the instant t = at the point xi is ip{xi), 
then the amplitude, iPt{x2), of finding it at X2 at time T will be 

/oo 
K{x2,T\xi,0)ijixi)dxi. (2.3) 
-00 

Here the propagator, K{x2,T\xi, 0), which describes the transition from xi-bol X2 in time T, is, 
says Feynman, a "sum" of the contributions (2.2) of all paths, 

P7. (2.4) 



K{x2,T\xi,0) = / exp 
Jr 



But all this intuitive, unless we say what "P7" means here. The definition of the integration 
measure is indeed the main stumbling block of the whole theory, and, despite many efforts, no 
fully satisfactory answer is available as yet. The miracle is that the integral can, in some cases, 
nevertheless evaluated [B]. Below we present one possible method we illustrate on the example 
of a one dimensional harmonic oscillator. 

Let us indeed consider an arbitrary path 7(i) in V that satisfies the given boundary condi- 
tions. Let us assume that V contains a unique classical path, 7(t), and let us decompose 7 into 
the sum of j{t) and a 

-fit) = 7(t) + r]{t), ahol r/(0) = r]{T) = 0, (2.5) 

since the end points are kept fixed. The Lagrange function of the oscillator is L = ^(x^ — w^x^). 
The action along 7(f) is therefore 

5(7)=/ \m{f - uj'^^'^)dt + m I {^^-uj^jrj)dt+ [ {ij^ - i^'^rj^)dt. (2.6) 
JO "'0 "'0 

For the oscillator, x + uj'^x = 0, so that x{t) = A sin ujt + B cos ujt., where the constants A and B 
are determined by the initial conditions, 

. X2- xi cos ujT 

B = xi and A = — , 

sm LoT 

provided sintoT ^ 0,i.e., if uoT ^ Ntt. Hence If T is not integer multiple of the half-period, 

T^Nx^, r=-, (2.7) 
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there exists a unique classical path that starts in xi and arrives, after time T, to X2-be. The 
action calculated along this path is 



2 sin ujT 



x\ + x\) cosujT — 2xiX2 



(2.8) 



Returning to (2.6), the first term is 5(7). Next, integration by parts of the middle term 
yields 



(7?) — uj'^-yri)dt 



+ uj'^^)r]dt = 0, 



since 7 satisfies the classical equation of motion x = —uj'^x. This term vanishes therefore, 
providing us with the propagator 



/r(x2, 0) = exp 



n 



S{l) 



X F(r), 



where the "reduced propagator", F{T), is a path integral taken over all variations. 



F{T) 



exp 



2h 



{rf^ - Lo'^rj^ dtj Vr]. 



(2.9) 



(2.10) 



But all this is still intuitive, as we still not say what "2?t/" actually means. To answer this 
question, Feynman expands the (periodic) variation into a Fourier series. 



= ^aksmlk-t] 
k=l 



f rj'dt = ^Y.-l f ^^dt=^E(^)'4. (2.11) 
Jo 2 ^ Jo 



Then he argues as follows. Any path is determined by its Fourier coefficients a^; instead of 
"integrating over all paths" let us integrate over all Fourier coefficients. 



where 



. . . / exp < i —al > X dai... dan 

"OO J-00 \,k=l J 

A. = m((^)2-^2) 



(2.12) 



(2.13) 



and where ^ denotes the Jacobian of the (linear) transformation V |Fourier coefficients |. 
Let us observe that J' is independent of the data (w, m, etc) of the oscillator and even of h. 
The classical Fresnel integrals can be evaluated. 



exp 



t—x 
2 



dx 



e'4. 



(2.14) 



Then the sum in the exponent in (2.12) becomes a product, 

/ n N 

where Cn is the product of various (divergent) factors. 
The product of eigenvectors can be split into two parts. 



"1/2 



(2.15) 



fe=l k=l k=l 



(2.16) 
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The Euler formula says now that 

oo / 

TT ( 1 _ 

A;%2 



(^-TJ^2]=— (2-17) 

k=l \ 



SO that 

V sm LoT 

where C denotes the product of all ^-independent factors. 

Let us now remember that for ^ we get a free particle, and therefore 



This yields C, and inserting the action (2.8) the propagator becomes, at last, 

K(x2,T|xi,0) = ( J""" ) ^ X exp I '"""^ [{xj + xl) cos loT - 2x1x2] ] . (2.19) 
V27rzn,smu;j / t27ismLL'i J 

3 The phase correction 

Let us observe that if the "sin" in the denominator vanishes, the whole expression becomes 
meaningless. How should we then continue after such a singular point. The answer comes from 



the evaluation of the Fresnel integrals : Euler's formula, (2.14), is only valid for A > 0. But 



Afc>0 ^ 0<T<k^, (3.1) 

where r = lir/uj is the period. Before reaching the first half period, < T < r/2, all factors 
under the square root are positive and Feynman's calculation is correct. After the first half 
period (but before a full periode), however, i. e. for r/2 < T < r, the first factor under the root 
is negative, while all the other factors remain positive. 

Ai < 0, Afe > 0, k> 2. 

As a result, the propagator gets multiplied by 

-^ = -i = e-'^. (3.2) 

Thus, the phase of the propagator (and therefore also of the wave function) jumps by (— 7r/2). 
Similarly, after N but before (A^ + 1), half-periods, i.e. for 

Ax ^ <T< (A^ + l) X ^, (3.3) 



the first A factors in (2.15) become negative. The phase jumps, therefore, by A^ x (— 7r/2). The 
correct result is hence 

(3.4) 

f imto n r, 1 ^iiL\r 

wrr-- 7F\ ^ [[xi + xi) cosljT - 2x1x2] > xe 2 '\ 

{ 2n\ smwi | J 

We must admit that our argument has been somewhat sloopy: each factor could be i 
or — i, and choosing the second one would change every phase jump from —tt/2 to tt/2. Worse: 
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every jump can be chosen independently ! Which square root of (—1) has to be chosen ? During 
the first half-period, the question is irrelevant, since it is merely a global phase. But repeating 
it times it is relevant ! 



The formula that generalizes (2.14) can be derived through analytic extension. For any real 
A/0, 



(^)VVt A>0 



exp 



■.A 2 

^-x 



dx 



A 

27r 



(3.5) 



|i/2e-i A<0 



Euler's formula, (2.17), only holds for x < vr, and otherwise it should be replaced by 

,.2 



n 

k=l 



1 



X 



smx\ 

X 



X > 0, 



(3.6) 



This confirms the validity of our previous argument : passing through every half-peridod 
contributes a new negative A, and this changes the phase by tt/2. 

What happens for T = Nxt/2 ? The propagator plainly diverges , since sin ujT sin Ntt = 
in the denominator. To derive the correct result, let us use the (semi-)group property of the 
time evolution. 



t+t' 



UtoUf 



N- 

2 



\2N 



According to (2.19) 



'ip{xi)dxi, 



(3.7) 



(3.8) 



4 ' ^ '2ttTi' 

which is, essentially, a Fourier-transformation. But the square of such a Fourier-transformation 
carries any function into itself up to a change of its argument. Hence 

V^^r(x2) = e-f^V((-l)^X2), (3.9) 



I.e. 



K{x2,T\xi,0) = exp 



2 



X 6{xi 



(3.10) 



Going to D dimensions, the phase will jump by D x 7r/2. What happens if the oscillator is 
not perfectly harmonic ? One has to study higher-order terms [8j. 
The result can be extended to several similar situations. 

• For a forced oscillator driven by a constant external force [6l[9l[T0], the previous calculation 
can be repeated word-by-word. The Lagrange function is 



L 



m 



[X — UJ X ) + JX 



(3.11) 



where / = const. For T ^ A;r/2 the Hamiltonian action is 



2 sin ujT 



{x\ + X2) cosojT — 2xiX2 



+2/ 



(1 — cosujT), , „r,2(l — cos LoT) — LoT sin LoT 

-{X1 + X2)- r- 



(3.12) 

cf. (2.8). The propagator is again (|3.4|) with the only change that the extra terms in (3.12) 



should be accounted for. In fact [10] . 

Kf{x2, T\xi, 0) = exp 



Kosc{x2 - x*,T\xi - x*,0) 



(3.13) 



where x* = f /muj^. This latter formula also holds for T = N x t/2. 
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4 Derivation from the wave function 



Before further investigating the various aspects and derivations of the Maslov correction from 
the Feynman integral viewpoint, we would like to show how it can be understood using our 
knowledge of the solution of the Schrodinger equation. The clue [10] is to write the propagator 
as 

oo 

K{xi,X2;T) = ^exp[-i{n+ ^)ujT]ipn{xi)ipnX2, (4.1) 

n=0 

where the tpn s are the normalized harmonic oscillator wavefunctions : 

where C = \/muj/fix and the Hn denote the Hermite polynomials. Their property important 
for us is their behavior under space reflection, il^ni—x) = (—1 )""(/'« (2^); implied by that of the 
Hermite polynomials. 

(4.1) can now be evaluated. Let us assume that we are not in caustic, t 7^ Ntt/uj, where N 
is some integer. Then write ujt = Ntt + ujt, < r < tt/uj. N is hence the integer part of wt/vr. 
Now from (4.1), we infer that 

00 

K = exp{—^iNTT) exp[— inA^TT — (i(n + |)u;r]-^„(xi)'0„(x2). (4-3) 

n=0 

00 

K(xi,X2;r) =exp(-|i7riV)^expH(n+|)u;t]V'„(xi)V„((-l)^X2). (4.4) 

n=0 

The Hermite polynomials admit a generating function, namely 

a(Ci,C2;.) = ^^exp y y (4.5) 

Q is analytic over the z plane with cuts going from —00 to —1 and from +1 to +00. For \z\ < 1 
the generating function admits the expansion 

00 n 

g{Cl,C2;z) = J2^Hn{Cl)Hn{C2). (4.6) 

The product of the Hermite polynomials is expressed therefore, as 



Then putting z = exp(— az), < a < vr, a lengthy calculation [TU] yields 



(4.7) 



-i7r/4 

g = exp[i(ia + Cl + (C2)2)]-— exp ( — — [C? + {Cifcosa- 2C1C2] ) • (4.8) 

2 sm a \2 sm a 

Collecting our formulae. 



K{xi,X2,T) = Jniu;/7rnexp{-lTrN[-]exp[-i{iu;T + Ct + Ci)]giC,{-irC2;e-'n. (4.9) 



vr 



Noting finally that sin(ti;r) = |sin(u;T)| = (— l)'^ sin(ti;t) and cos(a;r) = (— 1)'^ cos(wt) we get 
precisely the propagator ( 2.19| ), valid for times which are not integer multiples of half-periods. 
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Figure 1: Those paths which join xi to X2 in time T form an "infinite dimensional manifold" V. 
The the actual motion, 7(t), is a critical point of the hamiltonian action, viewed as a real- valued 
function on V. 



In a caustic i.e. if t = {Tr/uj)N, then, letting ujt Ntt in (4.4) yields [TO], using the 



completeness of the wave functions, the formula valid in the caustic point, i.e., 

lim K{xi,X2;t) = exp(-iiiV7r)(5(xi - (-1)^x2) (4.10) 

a;t— ^ A''7r 



cf. ( |3.10 ). 

Let us mention that the Maslov correction can also be derived from a semiclassical analysis 
of the wave function [2] . 

5 How minimal is the "minimal action" ? 

Let us now return to classical mechanics. The classical hamiltonian action is, as explained in 
Section [2] a real function, S{'y), defined on the "infinite dimensional manifold", V, of all paths 
which join xi to X2 in times T. A variation r] can be viewed in turn as a "tangent vector" to V 
in 7 [11], cf. Fig. 1. 5S, the first variation of the action, is the directional derivative of S in the 
"point" 7, 

5S,{v) = lim^^^^±^^'l^^. (5.1) 

s-^O s 

SS-y is, hence a one-form on V. Just like in finite dimensional calculus, if 7 is an extremal point 
of S, then the directional derivative must vanish in any direction, 

d/dL\ dL- 



6S^ 



This yields the classical (Euler-Lagrange) equations, ^(^ff ) ~ ~ 0' "^hose solution satisfying 
the boundary conditions, 7 is, by Hamilton's Principle, is the actual motion. 

But does such a solution always exist, and if it does, is it unique ? The answer is, generically, 
affirmative. Not always, however. This point can again illustrated by the oscillator. Let us 
indeed assume that the elapsed time is an integer multiple of the half-period, 

T = Nx^. (5.3) 

Then, independently of the initial velocity, all motions which start in xi arrive, at time T, to 
X2 = (-1)^X1. (cf. Fig. 2.). 
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1.5 




-1.5 I ' 1 1 ' ' 1 

0.2 0.4 0.6 0.8 1 

Figure 2: After a half period, all trajectories meat again, independently of the initial velocity, 
in the same point opposite to the one they started from. 



For this particular time, (i) either i.e. when X2 7^ (— there is no classical path at all 
with the required boundary condition, (ii) or i.e. for X2 = (— there are infinitely many 
of them. 

The solutions of the classical equation of motion behave, hence, as light rays that start 
from one focus of an elliptic mirror: after reflection, they are all collected into the same point, 
namely the other focus. (By analogy, such point are also called focal points also in variational 
mechanics) . 

Apart of this particular situation, there is a unique classical path, 7, in V. Let us assume 
that we are in this, generic, situation. Is the action minimal ? Let us emphasized that, just 



like in finite dimension, (5.2) is merely a necessary condition for having a minimum. It is not 
sufficient, though. In other words, 7 is a critical point of S, but not necessarily a minimal one. 

Let us mention, at this point, that the variational calculus is in fact local, i.e., it is a 
differential calculus in the neighborhood of a path 7. When calculating the first variation, 
we compare in fact the value of S calculated along paths which can be smoothly deformed to 
7. All such paths belong to a single path component of P. If the latter has more then one 
path components, as it happens in the, one has to perform the variational calculus in each 
path-connected sector separately. 

In technical terms, a "point" (i.e. a path in a manifold M) 7 belongs to the path-component 
of 7 if 7 can joined to the "point" 7 in P if there exists a continuous "path" of "points" 7^ such 
that 71 = 7, 70 = 7- But this means precisely that 7 can be deformed into 7 by a homotopy. 
The path components of V are therefore labeled with the homotopy classes of the underlying 
space, vro('P) = 7ri(M). 

An illustration is provided by the Aharonov-Bohm experiment [12], where M is the punctured 
plane R^\{0}, whose vri is Z. The (free) action has, hence, a minimum in two homotopy classes, 
namely in those which pass near the solenoid (assumed infinitely thin) on either side. In all other 
path components, i. e. for paths which turn around the origin, the action has no critical point. 

Returning to the critical points of the action, they behave just like those of a function defined 
over a finite (D) dimensional space : xq is a critical point of a function V{x) if 

6V{xo) = i.e. diV{xo) = V i. 

Whether a given critical point xq is a minimum, a maximum or a saddle point, can be determined 
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by looking at the matrix formed from second-order partial derivatives, 

Is the matrix of 5^V definite ? Being symmetric, will have D (real) eigenvalues, 

ea = Xae-a, a=l,...,D. 

Then xq is 

a minimum if > Va 

a maximum if Aa < Va (5-5) 

saddle point if Aa > A^ < for some a, b. 

In an analogous way, to see if a given classical path 7 makes the action minimal one or not, 
we have to calculate the second variation, This is a quadratic form on the "tangent space" 

of the "variations" . To be a minimum, the second variation has to be positive definite, 

d^Sj{r],r)) >0 Vr?. (5.6) 

If, however, 

6'^S^{ri,7]) <0 but 6'^S^{7j',r]') >0 (5.7) 

for suitable variations r] and rj' , then 7 is a saddle point. 
A variation 77 such that 

6^S^{r],r]) <0 (5.8) 

is a negative mode, and if 

ld^S{rj,r^')=0 Vry', (5.9) 

it is a zero-mode. In the direction of a zero mode the function is, in the first approximation, 
invariant. In this case, the nature of the critical point depends on the higher variations. 
Let us examine the second variation in some detail. 

where the integration is along the (assumed unique) classical path 7. After partial integration, 

where A is the operator of the second variation. \5'^S is hence positive definite if all eigenvalues 
of the quadratic form A are positive, A > 0, for all 

Ar? = A?7, 7/(0) = r/(r) = 0. (5.11) 

• In the simplest possible case of a one-dimensional massive particle in a potential, 

L= — -Vix) A = -7nU, + ^ . (5.12) 
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For the ID oscillator, e.g., 



A. 



(5.13) 



Hence, taking into account the boundary conditions. 



rj{t) = sin[(a;^ + 



a/2 



m 



t]. 



But TJu)^ ~^ m ~ ^'^^ k = 0, ±1, . . . due to periodicity, so that 



m 



ijj 



and r/fc(t) = sin(— t). 



(5.14) 



cf. (2.11) and (2.13). The integer iV > in (3.3) Nt/2 < T < {N + l)r/2, counts hence the 



negative eigenvalues. Thus, the oscillator trajectory is a minimum of the action during the first 
half period. For time beyond Nt/2, it becomes a saddle point with N negative modes ! 

What happens in the focal points ? For T = N x t/2, as seen above, either we have no 
classical path at all or we have infinitely many between our fixed endpoints. Let X2 = {—l)^xi 
and let us consider some 7. Then all the other classical paths can be viewed as a variation of 
7, labeled by a parameter s, 7 = 7^. The action calculated for any such classical path 7 = 7^ is 
the same. Varying s we get, hence, a "curve" in V inscribed onto the "level surface" S = const. 
The derivative w.r.t. s of S{js) is therefore zero, so that 



ds 



(5.15) 



is a null mode at 7. 

Minimum or not : what are the physical consequences ? At the purely classical level, nothing 



at all : (5.2 ) yields the correct equations of motion in all cases. In Quantum Mechanics, however. 



the consequence is precisely the Maslov phase jump, as we explain it below. 
Let us first consider another example. 

• A phase jump similar to the one found for the oscillator is observed for a charged particle 
moving perpendicularly to the induction lines of a constant magnetic field [6]. Classically, the 
particle rotates uniformly with Larmor frequency eB/m i.e. with period r = 2Tim/eB. The 
Lagrange function is 



^B = +y + 2uj{xy - yx) 



UJ 



eB 
2m 



(5.16) 



If T is not an integer multiple of a (full) period, T ^ N x t 



2Nnm 



, then there is a 

unique classical trajectory that links any two points (xi,yi) and (x2,2/2) in time T. The action 
calculated for it is 



Sb 



muj 



[{x2 - xif + (y2 - yif] coi{u:T) + 2{xiy2 - X2yi) 



(5.17) 



The matrix of the second variation is now 



Ab = m 



-2uj- 



dt^ ~~ dt 

d 

dt dt^ ^ 



2lo- 



(5.18) 
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Our task is to solve the eigenvalue problem 
mrix - 2mujr]y = -X-q^ 



ryx(o) = rjy{Q) = = r?,(r) = ny{T), (5.19) 



The solutions are readily derived by separating the equations (5.19) by applying a time- 
dependent rotation, 

r]x = cos u)t ^ + svciuot C, [ + (mu;^ + A)^ = 

r/y = — sin ^ + cos C [ mC, + {muP' + X)Q = ^ 

By periodicity, C oc sin {-s/uj"^ + \ t) = sin ^t, where k is some integer. The eigenvalues are, 
therefore, doubly degenerate, and are identical to those in the oscillator problem : 



m[(^f-u:^y A; = 0,1,.... (5.20) 



The reduced propagator is hence identical to that of a planar oscillator whose frequency is half 
of the Larmor value, uj = eB/2m. 
U Nt <T <{N + l)r, then 

K{x2,y2,T\xi,yiO) 



exp 



~2h 



2TTifi\ sin ujT\ 

(5.21) 

[(x2- xif + {y2- yif]cot{u)T) + 2[xiy2- X2yi) M-1)^, w = — . 

J I 2m 



cf. |13| [To] . After N full periods i.e. at T = ^^g'" [which corresponds to N half-oscillator- 
periods], all classical motions meet again in the point they started from. According to our 
previous results, taking into account the dimension of the problem, the propagator is again a 
Dirac-delta with a sign change : 

K{x2,y2,T\xi,yi,{)) = {-l)^ 5{xi - X2,yi - y2)- (5.22) 

It is worth mentioning that this calculation should turn out to be useful to explain the 
Sagnac effect [U]. This experiment, originally proposed and performed with light and later been 
repeated with massive particles, amounts to perform a two-slit type interference experiment, 
when the whole apparatus is fixed on a turntable. The clue is that the inertial force due to 
rotation behaves exactly as a fictions magnetic field, with twice the mass, 2m, replacing the 
electric charge. 

6 The Semiclassical approximation 

Returning to the general case, let us assume that V has a unique classical path, 7, and let us 
develop the action to second order : 

5(7) = 5(7) + <5^(r?) + 1525. rj) + ... (6.1) 

where the "dots, . . .", denote all higher-order variations. 

Taking the semiclassical approximation amounts of dropping all these terms. According to 
Hamilton's Principle 6Sy = 0. The semiclassical propagator is, hence. 



K{x2,T\xi,0) = exp 



xFiT), F{T) = y'exp|^,5?S(7?,r?)|p7?. (6.2) 
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cf. (|2^-( [2l0l ). 

The reduced propagator can be determined diagonalizing the second variation. Let us as- 
sume, for simphcity, that the system is 1-dimensional. The quadratic form A in the eigen- 
value equation (5.11) is a self-adjoint Sturm-Liouville operator on the space of all "tangent 
vectors" r]. The eigenvalues, A„, are therefore all real and form a complete orthonormal system 
w.r.t. the usual scalar product, {r]n,r]m) = J T]n'i]mdt = 6nm- Expanding the variation rj as 

r] = J^n^nVn yields (j],A7jj = J2k '^k^k- The reduced propagator is, therefore, once again (2.12), 
with the Xk denoting the eigenvectors of A in (5.11). Hence 

F{T) = C. 



2mfL 



(6.3) 



As explained before, the Jacobian is independent of the dynamics, so that (6.3) holds also for 
the free factor, 



\ 



n 



2i-K% 



free ' 
k 



(6.4) 



where F(r)^''^'= and yj""' are the free propagator and the eigenvalues, respectively. Dividing 
(|0| by K^ C drops out. 



F{T) = F{T)f'''^^ 
The semiclassical propagator is, therefore. 



K 



m 



2mhT 



1/2 



(6.5) 



since, in D dimensions, F{T)^^'^'^ = ( 



,D/2 



2'nihT> 



cf. (2.18). 



• For a ID oscillator, we recover the previous result : by (2.16), the quotient of the products 
of the eignevalues under the root is exactly the infinite product we determined using in the Euler 
formula, since A^'^'^^ = mvr^/c'^/r^. 



The Van Vleck matrix . 

Our result can also be presented using the Van Vleck matrix |15j . Let us indeed remember 
that, still assuming the uniqueness of the classical path, 7, between the to given points xi-et X2 
in time T, the action can be viewed as function of the end points. 

5(xi,X2) = 5(7) (6.6) 

is in fact Hamilton's Principal Function. The determinant of the DxD matrix 









dxidx2 




dx\dx2 



(6.7) 



is called the Van Vleck determinant. Then 



Theorem (|15]) : In D dimensions, the semiclassical propagator is 



K{x2,T\xi,0) 



1/2 



det 



dxidx2 



D/2 



exp 



-iDN- 



X exp 



-S(7) 



(6.^ 
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where xi = {x\) and X2 = {X2) are the initial and final point, respectively. 
• For the ID oscillator 



dxidx2 



muj 
sin ujt ' 



(6.9) 



so that (6.8) yields, once again, the previous oscillator propagator formula. 



Let us consider a charged particle in a constant magnetic field, and let us assume that 



T ^ T. Hamilton's Principal function is now (5.17), and the absolute value of the Van Vleck 
determinant reads 



Idet 



dxidx2 



rn^uj^ X det 




1 

cot LoT 



From (6.8) we get (5.21) once again. 



(6.10) 



7 Morse theory [ 1I6 



The methods of variational calculus allow us to further deepen our understanding. Let us chose 
a classical motion 7 [in D-dimensions] , and let us focus our attention to the second variation. 



As seen in (5.10 ) 



{r],Ar])dt, 



where A is the matrix of the second variation. A Jacobi field along 7 is a D-dimensional vector 
|(i) = {^i{t)) such that 

A^"'=0 i.e. A,j^j=0 \/i. (7.1) 

This second-order differential equation in 2D dimensions has 2D independent solutions. 
Generally, none of them belongs to the tangent space of V at 7, since ^(0) and/or ^(T) may not 
vanish. Two points p = 7(ii) and q = 7(ti) of a curve ^ , ^ < ti < t2 < T are called conjugate 
point, if there exists a Jacobi field such that 



(7.2) 



i.e., if it belongs to the tangent space. The number of such independent fields is called the 
multiplicity of the conjugate points p and q. 

The null space of the second variation is the set of all tangent vectors r] for which \5'^S{ri, r]') = 
for all 7]' . The dimension, u, of the null space is thus the number of independent Jacobi fields 
that also belong to the tangent space. The null-space is, hence, non-trivial, if the starting and 
the ending points of 7, p = xi = 7(0) and q = X2 = ^{T), are conjugate, v is finite and it is 
readily seen that v < D [16] . 

It is easy to see that 2D independent Jacobi fields can be constructed using classical paths. 
Let us indeed consider a one-paramater family of classical paths [i.e. solutions of the Euler- 
Lagrange equations] j3{u), < n < 1, such that /5(0) = 7. These paths are not required to have 
the same end points as 7. Then it is easy to see that 



dl3 



du 



u=0 



(7.3) 



is a Jacobi field along 7. The general solution of the Euler-Lagrange equations depends on 2D 
parameters. Deriving w.r.t. these parameters yields 2D independent Jacobi fields. 
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Let us now define the Morse index of 7 as tlie maximal dimension, /x, of those subspaces of 
the tangent space upon which the restriction of the second variation is negative definite. 

S^S^{ri,r]) <0. (7.4) 

Then we have : 



Theorem (Morse) [16] : Morse index is the number of those points of 7, each counted with its 
multiplicity, which are conjugate to xi = 7(0) . 

The essence of the proof is to show that the number of negative eigenvalues is precisely the 
Morse index fj, = fJ,(T). This is proved by showing that the Morse index, /u(t), calculated on the 
segment from t = to t of the curve is a monotonic function of t which, for sufficiently small t, 
is continuous from the left, /i(t — e) = fi{t) for e > sufficiently small. From the right, instead, 

fi{t + e) = fj.{t) + u, (7.5) 

where ly is, as before, the dimension of the null space of the second variation. The /x(t) as a 
function of time is hence, at first, constant, and then, upon crossing a conjugate point, it jumps 
by u. 

The Morse index is finite, so that our curve contains a finite number of points conjugate to 
the starting point. 

Let us again illustrate the general theory on example. 



For the ID oscillator and for A in (5.13), the solutions of the Jacobi equation (7.1) are 



a sin ujt + b cos cot, which vanish at t = and at t = T if & = and Tlo = Nit. A point q of the 
classical trajectory is conjugate to xi = 7(0) if 

q = qk = 7{k-), k = i,... 

Consistently with our general theorem, the null space of can be also obtained by considering 
the classical motion which starts from xi with velocity v, 

V dx(t^ 1 

— sin Lot + xi cos ujt =^ Cit) = — = — sinwt. (7.6) 

LO OV U) 

For t = we get manifestly a null-mode of the second variation. 

For a D-dimensional oscillator, each conjugate point contributes D null-modes. When T 
crosses 2N/t the dimension of the null space, i.e., the Morse index jumps by D. 

Let us mention that, on compact manifolds, the Morse index is related to the topology of 
the infinite dimensional manifold of paths [l6j . 



8 Conclusion 

Is the phase observable ? The naive answer to this question is negative: when the wave function 
is multiplied by any phase factor, the absolute value of the amplitude, and hence its square, 
the probability, are unchanged. This is only true for one wave function, though. If the sys- 
tem is decomposed into parts and the partial waves are multiplied by different phase factors, 
recombination will yield an interference. 

Similar phase jumps are found in molecular [T7], nuclear [18], and heavy ion [19] collisions. 
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After the original ideas reviewed here have been put forward, several developments took 
place. See, e.g. [201 [21]. Related questions are also discussed in some textbooks |22] . 

It is worth calling the attention to that all path integrals studied in this paper can be 
calculated [23, by transforming the system into the free form using a "non-relativistic conformal 
transformation" |24| |25] . For the harmonic oscillator, for example, every half-period can be 
conformally related to a full free motion. Then the Maslov correction can be recovered by 
fitting together the wave function at half-periods [25] . This latter paper also provides a complete 
catalog of all possibilities. 
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